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Abstract 

We show how it is possible to formulate Euclidean two-dimensional quantum gravity as the 
scaling limit of an ordinary statistical system by means of dynamical triangulations, which 
can be viewed as a discretization in the space of equivalence classes of metrics. Scaling 
relations exist and the critical exponents have simple geometric interpretations. Hartle- 
Hawkings wave functionals as well as reparametrization invariant correlation functions 
which depend on the geodesic distance can be calculated. The discretized approach makes 
sense even in higher dimensional space-time. Although analytic solutions are still missing 
in the higher dimensional case, numerical studies reveal an interesting structure and allow 
the identification of a fixed point where we can hope to define a genuine non-perturbative 
theory of four-dimensional quantum gravity. 
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1 Introduction 



In rf-dimensional quantum gravity the basic observables are transition amplitudes between 
rf-l-dimensional manifolds as well as correlation functions between "objects", e.g. scalar 
fields, separated by a geodesic distance r. We will discuss in detail why the last statement 
makes sense even if quantum gravity deals with dynamics of fluctuating metrics. Presently 
no theory of four-dimensional quantum gravity exists and it is not clear if one should 
attempt to formulate such a theory entirely in field theoretical terms or view it as being 
embedded in a larger theory (like string theory). As long as the question is not settled 
we should explore both possibilities. In this article we will review an approach based on 
conventional field theoretical methods, as they are known from the interplay between field 
theory and the theory of critical phenomena. The basic idea is to perform a discretization 
of the theory and in this way make it well defined. In the parameter space of the discretized 
theory we search for critical points where the "lattice spacing" can be taken to zero and 
contact made to continuum physics. The discretization is called "dynamical triangulation" 
and consists of a replacement of the functional integration over all equivalence classes of 
metrics with a summation over all abstract triangulations of the given manifold. The 
"lattice" spacing is the length of the links in the triangulations. It is of utmost importance 
that the discretization can be viewed as a discretization directly in the space of equivalence 
classes of metrics and not in parameter-space. 

We will show that it is possible to formulate Euclidean two-dimensional quantum 
gravity as the scaling limit of an ordinary statistical system by means of dynamical trian- 
gulations. Further, scaling relations can be derived and the critical exponents have simple 
geometric interpretations. Of course two-dimensional quantum gravity is much simpler 
than four-dimensional quantum gravity. Nevertheless the study of this theory is very use- 
ful since the observables defined above can be calculated analytically. In addition it is 
possible to test numerical methods, often very useful in the study of critical phenomena, 
and verify that they work well in two-dimensional quantum gravity. Dynamical triangula- 
tions can be used to define a theory of gravity in four dimensions too. In four dimensions 
we can not solve the model analytically, but it is possible to use numerical methods to get 
a rather complete picture of the phase diagram of the discretized theory. 



2 Definition of two-dimensional gravity 
2.1 Continuum notation 

The quantum partition function of two-dimensional quantum gravity can formally be 
written as 



Z(A,G)= J2 J . ZuA^l e~ s ^ A ^ [ V g <f>(Q e~ s ^. 



xeTop 



Vol(diff) 



where S e h[g; A,G] is the classical action 

S eh (g;A,G) = f d 2 ^ 



A 



4ttG 



(2.1) 



(2.2) 



while S m is the action for some matter fields coupled to gravity. In a full theory of 
gravity we might have to sum over topologies but it is not yet clarified how to do that 
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except in a perturbative expansion in the genus of the manifold. In the following we will 
therefore restrict the topology to be spherical, but allow for boundaries on the manifold, 
i.e. we consider spheres with a number of "holes". For fixed topology the two-dimensional 
curvature term will play no role since it is a topological invariant and we will ignore it, as 
well as the matter term, since we will concentrate on the geometric aspects of pure gravity. 

Let us write the formal expressions for the objects of interest in pure two-dimensional 
quantum gravity. The first are generalized Hartle-Hawking wave functionals, i.e. the 
amplitudes 

for n disconnected one-dimensional universes of prescribed length Li i = l,...,n. The 
second class of observables refers more directly to the metric structure of two-dimensional 
quantum gravity. We can ask about reparametrization invariant correlation functions. 
The simplest example is: 

A) = L VoTpf) e " %;A1 / / *ty/M)Wy/g(?) - R) (2-4) 

where d g denotes the geodesic distance with respect to the metric g a b(0 an d A is the 
cosmological constant. It is possible to define a whole set of such reparametrization invari- 
ant correlation functions of the geodesic distance by multiplying the densities \/g(^)d 2 ^ 
an( i Vdi^')^ 2 ^' m (2-4) by invariant tensors like 1Z(£) and 1Z(£'). Note that with the 
definition (2.4) geodesic distance becomes a meaningful concept even in quantum gravity, 
despite the fact that we integrate over all metrics. 

In certain situations it is convenient to consider universes with a fixed volume V . 
Corresponding to (2.1) we can write 

and in this class of universes the analogue of (2.4) is 

G ( R > r > = Lwdk) I J d2 t^) d2 t'^) WAte) - r)s(J d 2 tjm - v). 

(2.6) 

It follows directly from the definitions that 

poo rc+ioo A A 

Z(A) = dVe~ AV Z(V), Z(V)= —e AV Z(A). (2.7) 

JO Jc-ioo 27T2 

and similarly 

poo rc+ioo A A 

G(R;A)= dVe- AV G(R;V), G(R;V) = — e AV G(R;A). (2.8) 

JO Jc-ioo 27T2 

Note that the short distance behavior of G(R; V) is related to the (internal) Hausdorff 
dimension dh of the ensemble of manifolds defined by the partition function Z[A] since it is 
proportional to the average volume in a "spherical shell" a distance R from an arbitrarily 
chosen point £ as long as R <C V 1 ^ dh : 

G{R;V) ~ R dh ~ l for R<^V 1/dh . (2.9) 
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2.2 Discretized notation 



In two dimensions we have two local invariants, the area element d 2 ^^Jg and the local 
curvature 1Z(£). To discretize our manifold we use as building blocks equilateral triangles 
with side length a, i.e. area \^3a 2 /4. We can form an ensemble of piecewise linear mani- 
folds by gluing together the building blocks in all possible ways compatible with the given 
topology. In this way the order of a given vertex, i.e. the number of triangles which share 
the vertex, is not fixed. This degree of freedom matches perfectly with the reparametriza- 
tion invariant degree of freedom given by the local curvature 1Z(£). In fact Regge calculus 
is precisely a prescription which assigns curvature to piecewise linear manifolds in a way 
which agrees with the natural concept of parallel transportation on such manifolds. Ac- 
cording to Regge calculus [1] the curvature on two-dimensional piecewise linear manifolds 
lives on the vertices and to vertex i one assigns 

UidAi = 2(2tt - ]T A9 t (i)) = 2e t . (2.10) 

t3i 

The summation is over all triangles t which share vertex i, and A0t(i) is the angle of 
triangle t with vertex i. If the sum of the angles is 2ir, the neighborhood of vertex i is flat. 
The difference from 2ir is called the deficit angle. dAi stands for a local area element we 
can assign to the vertex, i.e. it is intended to be d 2 £y/g(£) in the neighborhood of vertex 
i. If At denotes the area of triangle t it is natural to define 

dAi = \^A t (2.11) 

6 t3i 

if we share the area of a triangle equally between its vertices. With these definitions we 
have for a closed piecewise linear manifold: 

^2 T^i dAi = 47rx, ^2 dAl = area ' {2.12) 

i i 

where \ denotes the Euler characteristic of the manifold. In our case, where all triangles 
are equilateral these formulas simplify a lot. If n 4 - is the order of vertex i and Nt the total 
number of triangles we have: 

2tt /3 
UidAi = —{6 - m), J2 dA < = ^ 2 N T . (2.13) 

i 

In the following it is convenient to rescale a such that J2i dAi = cl 2 Nt- 

The discretized version of the quantum partition function will be [2, 3, 4]: 

ZM = E 7^ e ~" NTa2 ( 2 - 14 ) 
TeT 

Here the summation is over a suitable class of (abstract) triangulations T and Cy denotes 
the symmetry factor of the triangulation T. This factor is present for closed surfaces for 
the same reason as the special symmetry factors for Feynman vacuum diagrams. Cy is 
equal to the order of the automorphism group of the graph T. In the following we will 
usually choose a = 1, dAi = unless explicitly mentioned otherwise. 
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Figure 1: A triangulation with a two-loop, the associated two vertices and the two triangles 
adjacent to one of the links in the two-loop. 



Until now we have not specified precisely which class of triangulations to use, but 
just stated that we glued the equilateral triangles together to form (closed) manifolds of 
spherical topology. If we blindly glue N triangles together along the links in order to 
form a closed manifold we can create somewhat pathological situations where there are 
closed two-loops (or even one-loops) on the surface even if the building blocks are triangles 
(see fig. 1). Strictly speaking such graphs cannot be identified with a simplicial manifold 
according to the standard mathematical definition since the simplical neighborhood of a 
vertex is not necessarily a disk. Whether we allow such two-loops (or one-loops) or not 
should not be important for a continuum limit since these are just structures at the cut- 
off scale. It has been verified a posteriori that it does not matter precisely which class of 
triangulations we use, as long as they allow an identification of the topology of the given 
piecewise linear manifold. We can write (2.14) as 



where 7jv denotes the triangulations made of N triangles. We denote Z(pi) the grand 
canonical partition function and Z(N) the canonical partition function since we can view 
H as the chemical potential for creating additional triangles. Z(N) has the interpretation 
as the number Af(N) of triangulations in 7jv- This number is exponentially bounded [5]: 



As a consequence Z(pi) is well defined for ^ > fj, c and divergent for ^ < /2 C . The point fj, c 
is the critical point and we should define the continuum limit as ^ — > /2 C . In this limit the 
large N part of the distribution will dominate in the sense that we can write: 




oo 
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While the exponent 7 has the interpretation as power correction to the exponential 
growing number of triangulations and thereby determines the leading singularity of parti- 
tion function, it also determines the part of the fractal structure related to the tendency 
to branch into so called baby universes. Let us first define a baby universe as part of the 
piecewise linear surface separated from the rest by a minimal "bottleneck". The precise 
nature of such a bottleneck depends on the class of triangulations we use. If we consider 
only regular triangulations such a bottleneck will be a loop of three links. If we cut our 
piecewise linear manifold along this loop we will separate it in two. Usually this separation 
will be trivial in the sense that the loop will be just the boundary of a triangle and all we 
have done is to separate this triangle from the rest of the triangulation. However, in case 
the loop is not the boundary of a triangle belonging to the triangulation cutting along the 
loop will produce a separation into two non-trivial parts of the triangulation. The smaller 
part is called a "minimal bottleneck baby universe", abbreviated "minbu" [6], the larger 
part the "parent". 

A moment of reflection will convince the reader that the average number of baby 
universes of volume V in the ensemble of triangulations of total volume N will be given 
by 

W(V)) N « -^VZ{V) (N - V)Z(N - V) (2.18) 

where 3! is the number of ways one can glue the two boundaries of the minbu and the 
parent together. The additional factors V and N — V reflect the fact that the minbu 
and the parent both have a minimal boundary i.e. for a generic large surface of volume 
V and no accidental symmetry factors there will be V such manifolds for each closed 
manifold since the boundary can be placed at any of the V triangles. If we assume that 
the canonical partition function is given by (2.16) we get: 

(N{V)) N ~ N [V(l - V/N)\<- 2 . (2.19) 

which has the obvious interpretation that the probability density for branching to a baby 
universe of volume V is V"'~ 2 (l - V/N) 1-2 . 

2.3 Scaling relations 

In addition to the entropy exponent 7 we can introduce the critical exponents v and r] 
known from the theory of critical phenomena: v is the exponent for the inverse mass or 
the correlation length and r] is the anomalous scaling exponent. Further, 7 can be given 
an interpretation as a susceptibility exponent. 

In order to define these exponents [7, 8, 26] let us first consider the discretized analogy 
of the continuum two-point function G(R;A) given by (2.4). For each surface, built of 
equilateral triangles, we have by Regge's prescription a metric. This means that we can 
define the geodesic distance between any two points on the piecewise linear manifold. 
Here we use a simplified definition instead. The geodesic distance between two triangles 
is defined by considering the path connecting centers of triangles. A geodesic path is 
one where the length is minimal, and the length is usually counted in units of the lattice 
spacing a. In the same way we can define the geodesic distance between two links as 
the shortest path along a sequence of neighboring triangles which connects the two links 
and we can define the geodesic distance between a link and a set of links as the shortest 
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geodesic distance to a member of the set. Let us now consider the sub-ensemble of surfaces 
where the geodesic distance between two links l\ and 1% is fixed to be r (which in units of 
the lattice spacing a is an integer). Denote this ensemble of surfaces by T(2, r). We define 
the two-point function by: 

oo 

G4r)=J2e-» N £ 1- (2-20) 

N=l TeT N (2,r,N) 

If we have two marked links the combinatorial factor Cy = 1. G^r) is the discretized 
version of the volume- volume correlator 2 G(R;A). 
We have the following theorem: 

Theorem: The two-point function falls off exponentially for r — > oo. 

lim ~logG>(r) = > Q 

r-»oo r 

In addition the mass m(pi) > for fj, > fj, c and m(pi) is a decreasing function of fj,. 

This trivial but important relation follows from the fact that 

G^ ri + r 2 )>G^ ri )G^r 2 ), (2.22) 

simply because each term on the rhs of eq. (2.22) can be given an interpretation as a 
term belonging to the lhs of eq. (2.22). This is illustrated in fig. 2. Eq. (2.22) shows that 
— logG |U (r) is sub-additive and this ensures the existence of the limit (2.21). In addition 
Tn(pi) > because G^r) is a decreasing function of r. Again this follows from general 
arguments which allow us to bound the number of triangulations with N triangles and 
two marked links separated by a distance r in terms of the number of triangulations with 
N triangles and two marked links separated by a distance r' < r. The same kind of 
arguments leads to the conclusion that m'(pi) > for fj, > fj, c , i.e. m(pi) is a decreasing 
function for fj, — > fj, c . 

Consequently we expect the following generic behavior of the two-point function: 

G^r) ~ e - m ^ r , r>l/m(/x); (2.23) 

G^r) ~ r l ~\ l<r<l/m(/x); (2.24) 

x(fJ-) = ^G>( r ) ~ (a* — Mc) -7 + less singular terms. (2.25) 

r 

We use the exponent r l ~' n for the short distance behavior in order to be in accordance 
with the general notation: 

GW-^ 1 ^, (2-26) 

where the first factor is due to the angular average over a spherical shell of radius r. Note 
that (2.25) gives an alternative interpretation of 7 as a susceptibility exponent. Recall that 



2 It is possible to associate an invariant volume 2a 2 /3 with each link since it has two adjacent triangles 
and each triangle contains three links. Alternatively we could have defined the two-point function using 
directly marked triangles instead of marked links. Since we are later going to consider loop-loop correlation 
functions the first definition is technically more convenient in two dimensions. 
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Figure 2: The inequality (2.22). Two triangulations with marked links separated by 
distances r\ and can be glued together to a triangulation where the marked links have 
a distance r\ + but the same number of triangles by cutting open a marked link in each 
of the triangulations to a 2-loop boundary and gluing together the two boundaries. 



for a spin system the susceptibility is the second derivative of the free energy with respect 
to the external magnetic field and it has the alternative interpretation as the integral of 
the spin-spin correlation function. The same is true here: By (2.15) and (2.17) it follows 
that the second derivative of Z(pi) with respect to ^ has a leading singularity (pi — // c )~ 7 , 
while eq. (2.25) provides us with the alternative interpretation as the integral of the 
volume-volume correlator. x(m) ~ ^"(i 1 ) since both can be viewed as the summation over 
all triangulations with two marked triangles or links. 

Although we have proven above that the mass decreases as ^ — > /2 C the general argu- 
ments do not allow us to conclude that m(ji) scales to zero at the critical point /2 C . Let us 
add this as an assumption (we will later verify that it is true in two-dimensional quantum 
gravity) : 

m{n) ~ (fi-fi c Y, (2.27) 

and explore the consequences. In order to make the scaling relations more transparent let 
us for a moment reintroduce the dimensionful lattice spacing a and define the renormal- 
ized cosmological constant A in terms of the bare cosmological constant fj, by an additive 
renormalization: 

H = i_i c + Aa 2 , i.e. a(/i) ~ (fj, - n c ) 1 ^- (2.28) 

If we want a scaling limit where the concept of a mass survives we have to introduce a 
scaling like: 

m(ji) = M [a(ji)] 2lJ , i.e M = cA" , (2.29) 

R = r[a(ii)] 2l/ 1 (2.30) 
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where c is a constant of order one. In this way the scaling limit is one where the "bare" mass 
m(n) scales to zero as the lattice spacing vanishes, while the "renormalized" mass M is kept 
fixed. Similar remarks apply for the "bare" geodesic distance r and the "renormalized" 
geodesic distance, only will r go to infinity for a(pi) — > 0. The long distance behavior of 
the two-point function will be 

G^r) ~ e~ MR for fi fi c , m(/z) r > 1. (2.31) 

The relation between the continuum propagator G(R; A) and G^r) is unambiguously fixed 
by (2.23)-(2.25): 

G(R;A)= lim [a{n) 2v ] l -' n G^r), m(fi)r = MR. (2.32) 

In addition the only possibility for this limit to exist is that the long distance behavior of 
Gfj,(r) is given by 

G,(r) = /(m(/z)r) K^T" 1 e~ m ^ r , (2-33) 

where f(x) (which we ignored in (2.23)) is subleading compared to the exponential func- 
tion, i.e. 

G(R;A) ~ R 1 - 11 , cA^Cl, G(R; A) ~ f(ck v R) A^ -1 ) e~ cAl,R , ck v R > 1. (2.34) 

We shall verify that we indeed have the correct asymptotic behavior (2.32) (with f(x) = 
const, and v = 1/4) in pure two-dimensional quantum gravity. 

In the scaling limit (2.31) the exponent v has the simple geometric interpretation as 
the inverse of the internal Hausdorff dimension of the ensemble of piecewise linear surfaces 
T(2,r). Let (•),, denote the average with respect to the ensemble T(2,r). We define the 
internal Hausdorff dimension dh by 

(N) r ~ r dh for v — y oo, m(/i) r = const. (2.35) 

From (2.23) we have 

{N)r = J^G,{r) = d_rn^^ 

and using (2.27) we get 

(N) r ~ r 1/l/ , i.e. d h = l/v. (2.37) 
Another consequence of the above scaling is that we can rewrite (2.36) as: 

(N) R ~^-. (2.38) 

This relation shows how the number of triangles diverges as fj, — > fj, c , (as we want in order 
to connect to continuum physics) and in terms of the dimensionful volume V = Na 2 we 
have: 

(V)r ~ cA»R » 1. (2.39) 

This relation simply tells us that the typical universe which for a given cosmological 
constant A has points separated by a "continuum" geodesic distance R which is much 
larger that l/A 1 ~ 1 / dh will be a long tube. 
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It is interesting to give a direct physical interpretation of the short distance behavior 
of Gfj,(r) as defined by (2.24). In order to do so let us change from the grand canonical 
ensemble given by (2.20) to the canonical ensemble defined by T(2, r, N), the class of 
triangulations with N triangles where two links (or triangles) are marked and separated 
by a distance r. On this ensemble we can define the discretized analogue of G(R, V): 

G(r,N)= J2 !. ( 2 - 4 °) 

TeT{2,r,N) 

For r = we have the following N dependence (for the r dependence see (2.45) below) 

G(0,N) ~ N"<- 2 e^ N . (2.41) 

The partition functionary) ~ N 1 ~ 3 e flcN and the one-point function for large N is pro- 
portional to NZ(N) since it counts the triangulations with one marked triangle (or link or 
vertex depending on the precise definition). For r = (or just small) there is essentially 
no difference between the one-point function and G(0,N). 

G(r,N) is related to G^r) by a (discrete) Laplace transform: 

G^r) = J2G(r,N)e-» N . (2.42) 

N 

The long distance behavior of G(r, N) is determined by the long distance behavior of G^(r). 
Close to the scaling limit it follows by direct calculation (e.g. a saddle point calculation 3 ) 
that 

G^r) ~ e - r ^-^ 1/dh 



G(r,N) ~ e -c(r^/N)^ e , cN for r d h>Ni (2 _ 43) 

where c= (d h - I) / d d h hl(dh ~ l) . 

On the other hand the short distance behavior of G |U (r) is determined by the short 
distance behavior of G(r, N) which is simple. Eqs. (2.35) and (2.37) define the concept 
of Hausdorff dimension in the grand canonical ensemble. A definition in the canonical 
ensemble would be: Take N 1 ^ dh ^> r and simply count the volume (here number of 
triangles) of a "spherical shell" of thickness 1 and radius r from a marked triangle, sum 
over all triangulations Tjy with one marked triangle and divide by the total number of 
such triangulations. Call this number (n(r)) N . The Hausdorff dimension is then defined 

by 

(n(r)) N ~ r dh ~ l for 1< r < N 1/dh . (2.44) 

It follows from the definitions that we can write 

G(r,N) . 
(n(r)) N ~ i.e 

G(r,N) ~ r dh-i N -y-2 e ncN for 1 < r < N i/d h _ 



3 In addition to the exponentially decaying part of G(r, N) there is also a power correction coming from 
the quadratic integration in the saddle point approximation. We shall not consider the explicit form of the 
power correction here. 
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We can finally calculate the short distance behavior of G^r) from eq. (2.42). For fj, — > fj, c , 
we get: 

G,(r) ~ r^" 1 N^- 2 e-< rd ^ N )^ ~ r^" 1 (2.46) 
N 

This shows that 

r} = 2-jd h , i.e. 7 = 1/(2-77), (2.47) 

which is Fishers scaling relation, which is valid also for quantum gravity. Of course the 
relation could be derived directly from (2.23)-(2.25), but the above arguments highlight 
that the anomalous scaling dimension tj is a function of the two kinds of fractal structures 
we can define on the ensemble of piecewise linear manifolds: the Hausdorff dimension 
and the baby universe proliferation probability. In addition the arguments show that the 
canonical and grand canonical definitions of Hausdorff dimension agree. 

A most important remark is that the definitions and the scaling relations above gen- 
eralize to higher dimensional quantum gravity. 



2.4 Branched polymers 

The model of branched polymers (BP) provides us with a simple, but non-trivial example 
of the above scenario [10] and will play an important role in the following. In a certain 
way it can be viewed as the lowest dimensional fractal structure and it will appear as the 
limiting case of higher dimensional gravity theories. 

Let us define branched polymers as the sum over all tree graphs (no loops in the 
graphs) with certain weights given to the graphs according to the following definition of 
the partition function: 

Z(») = '£ 7 ?-p(BP)e-»\ BP \, (2.48) 
bp Cbp 

where \BP\ is the number of links in a BP and ^ is a chemical potential for the number 
of links, while 

p(bp) = n ( 2 - 49 ) 

ieBP 

where i denotes a vertex, n 4 - the number of links joining at vertex i and /(n 4 -) is non- 
negative. f(ni) can be viewed as the unnormalized branching weight for one link branching 
into ni — 1 links at vertex i. Finally C'bp is a symmetry factor such that rooted branched 
polymers, i.e. polymers with the first link marked, is counted only once. 

This model can be solved [10, 11]. It has a critical point fj, c (depending on /) and close 
to the critical point we have: 

Z'V)-^-^)- 1 / 2 , (2.50) 

i.e. 7=1/2 for branched polymers. On the branched polymers we define the "geodesic 
distance" between two vertices as the shortest link path, which is unique since we consider 
tree-graphs. The graphic representation of the two-point function is shown in fig. 3. Had 
it not been for the ability to branch, the two-point function would simply be 

G,(r) = e~» r . (2.51) 
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Figure 3: The graphical representation of the two-point function for branched polymers. 
The dashed line represents the unique shortest path between the two marked vertices. The 
"blobs" represent the contribution from all rooted polymers branching out from a vertex. 



However, the insertion of one-point functions at any vertex leads to a non-analytic coupling 
constant renormalization and the result is changed to [10] 

G^r) = const. e^vT^ for ^ _ ^ 0) ( 2 .52) 

where k is a positive constant depending on /. We can now find G(r, N) by an inverse 
Laplace transform: 

G(r, N) = const. N~ 3 l 2 r e -^ 2 /4iV (2.53) 

We confirm from this explicit expression that the (internal) Hausdorff dimension of branched 
polymers is 2 (like a smooth surface !) and that 7=1/2 since the prefactor of G(r, N) for 
small r should be 7V 7_2 r <iff_1 . 

2.5 Perspectives 

We have derived general scaling relations for quantum gravity. In the following we will 
solve the theory explicitly in two dimensions. It is remarkable that both the Hartle- 
Hawkings wave-functionals and the volume-volume correlation functions can be found by 
purely combinatorial methods at the discretized level, after which the scaling limit can be 
taken unambiguously. 

3 Matrix models 

In this section we will discuss how to construct Hartle-Hawkings wave functionals in the 
discretized approach. The problem was solved already in 1963 by the mathematician 
Tutte, who found the generating functionals for the number of planar triangulations with 
boundaries. We shall rederive here the results using the so called matrix models technique, 
following the approach in [12, 13]. 
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3.1 The Hartle-Hawking wave functionals 

In order to count the number of triangulations we represent the triangles by means of 
Hermitian matrices: Label the vertices of the i th triangle by abstract indices a>i,/3i,yi and 
attach a Hermitian matrix 4> at f3 t to the oriented link from oi{ to /3;. In this way we can 
attach the scalar quantity 

(^a.p.^pn.^.a, = Tr^ 3 (3.1) 
to each of the K triangles. The Gaussian integral 

iTr^ W1 T ,3 



d4>e ~ K\ U Tr W (3 ' 2) 
where 

d(f> = Y[ dMt4> a p Y[ dQ(f> a /3 (3.3) 

a<f3 a<f3 

can be performed by doing all possible Wick contractions of (^-fields. This corresponds to 
performing all possible gluings of surfaces of K triangles, the reason being that each Wick 
contraction in the Gaussian integral glues together two links: 

((kaprfajp,) = J d(f> e-^^ 2 <j) atl3t <j) ajl3j = 8 atl3j 8,3 taj . (3.4) 

This is illustrated in fig. 4. After all Wick contractions are performed in eq. (3.2) the 
K triangles have been glued together in all possible ways. The surfaces created in this 
way will consist of disconnected parts, but we get the connected graphs by taking the 
logarithm of all graphs. Furthermore we can calculate the contribution from a particular 
graph constituting a closed surface: we pick up a factor N, N being the number of indices, 
whenever a vertex becomes an internal vertex in the process of gluing together links by 
Wick contractions. This means that we get a total factor N v , where V is the number of 
vertices. If we make the substitution 

Tr^ 3 ^-^=Tr^> 3 (3.5) 



it is seen that the factors of N for K triangles combine to N v K l 2 = 7V X , since the Euler 
characteristic for a triangulation of K triangles, L links and V vertices is 

X = V-L + K = V- K/2. (3.6) 
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By multiplying each triangle Tvcf) 3 by 1/y/N the weight of a given triangulation only 
depends on its topology. In addition the sum over all triangulations exponentiates. Col- 
lecting this information we can write: 

where Z(fj,, G) is partition function defined in terms of dynamical triangulations and 

Z(g,N) = jd<f> exp (-^Tr^> 2 + ^Tr/) (3.8) 

provided we make the identification: 

i = log7V, fi =- log g. (3.9) 

Eq. (3.8) allows a 1/N 2 expansion and this expansion is an expansion in topology. Here we 
will only be interested in spherical topology, i.e. the leading term in the 1/N 2 expansion. 

For the purpose of a general (perturbative) analysis of the matrix integral (3.8) it is 
convenient to consider the generalization to an arbitrary set of coupling constants {</;}: 

Z( gi ,g 2 ,...) = J d^e~ NTrV ^ (3.10) 

where 

CO 

9 



v ({9i}) = E t*"- ( 3 - n ) 



1 n 

n=l 



In eq. (3.11) we have of convenience scaled <f> — > yN<f>. In this way the topological nature 
of the expansion is still preserved: All two-dimensional complexes of Euler characteristic 
X will have a factor N x associated with them. The interpretation of (3.10) is intended to 
be as before: we have in mind a Gaussian integral around which we expand, i.e. g 2 > 
and g n < with the sign convention used in (3.11). The convenience of considering an 
arbitrary potential is that the general coupling constants g n act as sources for terms like 
Trcj) n , and by differentiating Z with respect to g n we can calculate expectation values of 
"observables" like Tr^> ra . (Ti(f> n /N) has the following obvious interpretation: It represents 
the summation over all "surfaces" which have a ra-sided polygon with one marked link as 
boundary. This follows from the gluing procedure realized by Wick contractions of the 
Gaussian integrals 4 . In a similar way 

^<Tr<TTr<n - ^<Tr <T>(Tr <T> (3-12) 

will represent the sum over all connected two-dimensional complexes which connect one 
boundary consisting of n links with another boundary consisting of m links, i.e. precisely 
the discretized version of the Hartle- Hawking wave Junctionals. Let us define the generating 



4 In case we want unmarked links we should multiply (Tr <f) n /N) with an additional symmetry factor 
\jn. 
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functional for connected loop correlators. The expectation value of an arbitrary observable 
is defined by 

(Q(cf>)) = ±J # e - NW WQ((i). (3.13) 

The generating function for s-loop correlators, which we will also, somewhat inaccurate, 
denote the s-loop correlator, is defined by 

(3.14) 

where conn refers to the connected part as defined by (3.12), or its generalization to more 
correlators. If we introduce the so-called loop insertion operator by 

dV(z) ' = ~ "S^ 1 dfk (3 ' 15) 
it follows from the definition (3.14) that 

'"<*•••••'•> = ^wh)---^k) w{li) (3 - 16) 

and this equation shows that if the one-loop operator is known for an arbitrary potential, 
all multi-loop correlators can be calculated. 

The one-loop correlator is related to the density p(X) of eigenvalues defined by the 
matrix integral as follows : 

N 

p(A) = <5>(A-A,-)> (3.17) 

8 = 1 

where A;, i = 1, . . . , N denote the N eigenvalues of the matrix (f>. With this definition we 
have 

W(z) = r dX (3.18) 

J-oo Z — A 

For N — y oo there exist, as we shall see, consistent solutions where the support of p is 
confined to a finite interval [co, c\] on the real axis. In this case will be an analytic 

function in the complex plane, except for a cut at the support of p and it follows from 
Schwartz's reflection principle that 

2irip{X) = lim W(X + ie) - W(X - ie) (3.19) 

e— >0 

3.2 The loop equations 

Let us explore the invariance of the matrix integral (3.10) under field redefinitions of the 
type: 

<^^ + £]T^ r = ^ + £ -. (3.20) 

k^o p p ~ $ 

This kind of field redefinitions only make sense if p is chosen on the real axis outside the 
support of the eigenvalues of (f>. As mentioned above we will verify that this scenario is 
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realized for N — > oo. Under the transformation (3.20) the measure and the action change 
as 

d<f> d<f> ( 1 + eTr — !— ■ Tr —^—] (3.21) 

V p-<p p- <p) 



TtV(4>) Tr V(4>) + eTr ^—L^y'^ • ( 3 - 22 ) 



The integral (3.10) will be invariant under a redefinition of the integration variables by 
eq. (3.20) and the change of measure and action has to cancel to first order in e. By use 
of eqs. (3.21) and (3.22) this leads to the following equation: 

^(Tr-l^-^VTr (-L-V'tf))} e~ N ^ = 0. (3.23) 

The first term in this equation is by definition 

N 2 W{p)W{p) + W(p,p). (3.24) 

The second term in eq. (3.23) can be written as an integral over the the one-loop correlator 
by means of the density of eigenvalues and we can finally write (3.23) in the standard form, 
known as the loop equation [14]: 

where z is outside the interval [co, c\] on the real axis and the contour C encloses the 
cut, but not z (see fig. 5). In eq. (3.25) the 1/N 2 expansion is manifest and since we are 
interested in the leading term (spherical topology) we ignore the last term on the rhs of 
(3.25) and denote the corresponding solution Wo. 

In case of a Gaussian potential V(z) = z 2 /2 the eigenvalue density p(X) is given by 
famous Wiener's semicircle law: 

MA) = ^(2-A)(2 + A) i.e. W (z) = \ [z - y/(z - 2) (z + 2)) , (3.26) 

where the last equation follows from (3.18). 

For a general potential we can find a solution which has essentially the same structure 
as for the quadratic potential and where the eigenvalue density is given by [15]: 



/ o(A)~M(A)^/(c 1 - A)(A-c ), (3.27) 
where M(A) is a polynomial of order n — 2 if V(X) is a polynomial of order n. We write: 

Wo(z) = \ (v'{z) - M{z)^{z- Cl ){z-c )^ (3.28) 

and the requirement that Wq{z) = 0(1/ z) uniquely determines M(z) since (3.28) allows 
us to write 

[ du M(u) f du V'(lo) 1 
M(z)=<t — : —^-L=<f, — ; v ; w 3.29 

JCoo 2m z -lo Jcoo 2m z - lo ^(uj - ci) (u - c ) 
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Figure 5: The integration contour C and the cut from cq to c\. When deforming the 
contour to infinity we get two contributions: one from the circle and one from the 
circle C z around the pole z. 



as the part of the integral which involves Wq{z) vanishes. By expanding the last integrand 
in powers of — c\) we can write 

M(z) = Y J M k (z-c 1 ) k -\ M k [c ,c 1 ;g i ]= f — ^ (3.30) 

It is important to note that M k = for k > n where n is the order of the polynomial 
V(z). It follows directly from the integral which defines M k . 

Wo(z) can now be expressed in terms of the "moments" M k [12, 13]: 

Wo(z) = l - (v'(z) - y/(z - Cl )(z - c ) M k [c 0l Cl - l9t ](z - c ) fc -^ . (3.31) 

The endpoints x and y of the cut are determined selfconsistently by the following boundary 
conditions 

M_i[c ,ci;</ t -] = 2, M o [co,ci;</,-] = 0, (3.32) 

These two equations follow by contraction of the contour in the last integral in (3.29) to 
C. 

3.3 Complete solution Wo(zi, z s ) 

In principle the solution at spherical level is given by (3.31)-(3.32). These equations define 
Wo(z) and we can apply the loop inserting operator to obtain any multi-loop correlator. 
Quite surprisingly one can obtain an explicit formula if we change the matrix model 
slightly. Instead of Hermitian matrices we use general complex matrices and consider a 
potential: 

OO 

n 



V(cpU) = Y.-^(<pU) n . (3-33) 



n=l 
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The one-loop correlator is given by: 

n 

This model has a surface representation [16]. Tr ($<f>) n represents a 2ra-gon where the 
boundary links have alternating black and white colors, corresponding the <j> and ft. The 
effect of Gaussian integration with respect to the complex matrices is to glue together such 
"checker-board" polygons just as Hermitian matrices glued together ordinary polygons. 
The only additional rule is that only "white" and "black" links can be glued together since 
/ while ((f) 2 ) = and = 0. Such short distance differences in gluing should 

be unimportant in the continuum limit. 

We can write down the loop equations for this model [12]. They are similar to the 
ones already considered, except that we now have a symmetry with respect to (f> — > —cj), 
i.e. Co = — c\ and the formulas above simplify: 

M(z) = Y,M k [c ug% ](z 2 -c\) k -\ M k [c ug% ] = f £ 2 "_ j" j +1/2 , (3.36) 

k=l ^ 1' 

and the position c\ of the cut is determined by 

M [ Cl , </,-] = 2. (3.37) 

While these formulas look rather complicated it is a pleasant surprise that things sim- 
plify a lot for the higher loop correlators. In order to apply the loop insertion operator we 
note that 

d _ k d d c\ d 

dV(z) = " Y = 8V(z) ~ Mi {z 2 - Cf) 3 / 2 dc 2 ' (3 ' 3 ^ 

where 

d _ v k d dV'juj) 2ujz 

dV{z)~ Y z 2 ^ 1 dgk dV(z) z 2 -uj 2 ' { ' 

It is now straightforward, though tedious, to apply the loop insertion operator. We find 
[12]: 

1 



W (z 1 ,z 2 ) 



4(z 2 - z 2 ) 2 



2 \ Z \ C l , 2 Z 2 C l o 



(3.40) 



Wo(z 1 ,z 2 ,z 3 ) = -A- 1 (3.41) 

16Mx ^ {zl _ c r ){zl _ cl){zl _ cl) 

*> = G^) 3 2%ui n 2(4-4)*/* ' (3 - 42) 

Note that the above formulas are valid for any potential V . All dependence on the coupling 
constants is hidden in M\ and c\. 
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3.4 The scaling limit 

In the case of the simplest potential (Tr^> 3 for the Hermitian matrix model, Tr ((^(j)) 2 for 
the complex matrix model), we have one independent coupling constant g if we fix the 
coupling constant in front of the Gaussian term. Eq. (3.9) gives the relation between the 
bare cosmological coupling constant p, and the coupling constant g of the matrix models. 
We have seen that there is a critical p c such that the continuum limit should be taken for 
p —7- p c . Corresponding to p c there will be a g c . If we introduce the lattice spacing a, 
the area of our universe for a given triangulation is given by Nxa 2 , Nt being the number 
of triangles in T and as already discussed it is natural to introduce the renormalized 
cosmological constant A by 

p-p c = Aa 2 , i.e. g c - g ~ g c Aa 2 . (3.43) 

Let us for simplicity consider the complex matrix model. The density of eigenvalues 
p(X) is given by (3.27) and (3.36): 

p(X) ~ M{X)^c 2 - A 2 , M(X) = Mi + M 2 (c 2 - A 2 ). (3.44) 

The only non-analytic behavior is associated with the endpoints of the distribution where 
| A | — >- ci . If Mi[<7] is positive the behavior will be identical to that of the Gaussian model. 
We conclude that the critical point g c is determined by: 

M 1 [c 1 (g c ),g c ] = 0. (3.45) 

A glance on Wo(zi, ...,z s ) corroborates this observation since it will be singular precisely 
when M\ = 0. Let us move slightly away from the critical point by scaling g according to 
(3.43): 

g = g c (l - ha 2 ) (3.46) 

Corresponding to this change there will be a change c\' — > c\ ' —5c\. It can be calculated 
directly from (3.37) by expanding Mo[ci,</] and we get: 

(5c 2 ) 2 = Aa 2 + 0(a 3 ) i.e. c\ = c ( f - aVX + 0(a 2 ), (3.47) 

where we have absorbed a constant k 2 = 16/3M2(g c ) m the definition of the cosmological 
constant A. In addition Mi[ci,gi\ will now be different from zero: 

Ml [c 1 ,g t ] = ^^5c 2 + 0(a 2 ). (3.48) 

From the explicit formulas for the multi-loop correlators it follows that the complex vari- 
able z appears in the combination z 2 — x 2 and it is natural to introduce a scaling of 
z 2 : 

z 2 = z 2 + aZ, z c = c[ c) (3.49) 
With this notation we can write in the scaling limit : 

M(z) = Mi + M 2 (z 2 - x 2 ) = aM 2 (g c )(Z - ^VX) + 0(a 2 ) (3.50) 
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Z serves the same role in the scaling limit as z at the discretized level: knowing Wo(zi, ...z s ) 
allows us to reconstruct the multi-loop correlators consisting of discretized boundary loops 
of length rii, n s by multiple contour integration. In the scaling limit the physical length 
of these loops will be Li = n 4 -a, i.e. they will go to zero. If we want genuine macroscopic 
loops we have to scale n 4 - to oo at the same time as a — > such that Li is constant. 
By substituting (3.49) and (3.50) in the expressions for Wo(zi, ...z s ) we get an expres- 
sion Wo(Zi, Z s ; A) and we can reconstruct the corresponding multi-loop amplitude 
Wo(Li, L s ; A) by an inverse Laplace transform. The contour integration is changed 
into an inverse Laplace transform since the cut [— Co, c\] (or [co, c\] for the Hermitian ma- 
trix model) by the substitution (3.49) and (3.50) is changed into a cut [— oo,a/A]. The 
contour integration around the cut can now be deformed: 



Jc > r=i 



i z 2n 'W (z 1 ,...,z s )^ (3.51) 

c+«oo s 



„2(n 1 + -+n s ) 



/ ]JdZ t e L ' z ' W (Z ll ... 1 Z s - 1 A) 1 c>-VX 

J c—ioo 



(3.52) 

since we have 

z 2n = z 2 c n {l + aZ) L l a ~ z 2 c n e ZL 1 L = n/a. (3.53) 

The highly divergent factor z 2 ^ Ll ^ ' rLs )l a j s a W ave function renormalization of the macro- 
scopic boundaries. It is to be expected since it is possible in the continuum to add a term 



S M = A / dl = \ d£ </?(£) (3.54) 

JdM J 

to the action. This is just the induced one-dimensional gravity on the boundary of the 
manifold. Since A has the dimension of mass we expect in the discretized version that 
the bare coupling constant will undergo an additive renormalization (like the cosmological 
constant itself), i.e. we will have to cancel a term like z 2n . 

From eq. (3.42) we immediately get the generating functional for macroscopic multi- 
loop amplitudes [12]: 

W (zi,...,z s ) ~ a 5 - 7s ' 2 W(Z u ...,Z s ;A), (3.55) 

d s ~ 3 s 1 

W(Z 1 ,...,Z S ;A) = — — r TT 7=—r, s > 3- (3.56) 

The expressions for W(Z,A) and Z(Zi, Z^\ A) are slightly more complicated since Wq{z) 
contains a non-universal part (V'(z)/2). and we will not give them here. 
We can calculate the inverse Laplace transform (3.51) 

W{L u ...,L s ;k) = 

rc+ioo s r ]s-3 p -y/K(L 1 + -+L s ) 

/ Y[dZie ZiLi W(Z u ...,Z s ;A) 

J c—ioo ■ i 



i=l 



dA s ~ 3 yJL x ■■■L s 



(3.57) 



These are the Hartle-Hawkings wave Junctionals (2.3) of two-dimensional quantum 
gravity. 
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4 The two-point function 

Let us now turn to the problem of calculating the two-point function of two-dimensional 
quantum gravity as a function of the geodesic distance. 

4.1 Formulation of the combinatorial problem 

Let F(x;g) denote the generating functional for the number Af n ,i of triangulations of the 
disk with n triangles and one boundary of length /. It is given by F(x;g) = Wo(l/x; g)/x, 
where Wo(z; g) is the one-loop function for a cubic matrix model calculated in the last 
section 5 : 

oo 

F(x;g) = J2^n,W n x l = J2 X ' F (^S) 

n,l 1=0 

= K^ _ i) +/(s;5) ' 

where c 4 - are functions of g determined by the requirement 6 F(x;g) ~ 0(1) for x — > 0. 
We have F{l;g) = (Tr 4> l /N) in the matrix model notation of the last section and it 
has the interpretation as the generating functional for all triangulations with a boundary 
consisting of / links of which one is marked. For obvious reasons we denote F(x; g) the disk 
amplitude. As usual we have g = e _A1 where fj, is the cosmological constant. Corresponding 
to the critical value fj, c we have a critical value g c as explained in the last section and for 
fj, close to fj, c we have 

A/i = /i - /i c = Ag/g c . (4.3) 

One has cq < < c\ < ci as long as g < g c . The only thing we need to know here is that 
for g = g c we have c\ = ci (which we denote l/x c , while for g close to g c , i.e. fj, close to 
fj, c we have (as follows directly from (3.47)-(3.50) 

c 2 (/x) = l/x c + -yfAjl+0(An) (4.4) 

Cl (/i) = l/x c -ay/Ajl+0(An) (4.5) 
c (/x) = c (mc) + O(Afi) (4.6) 

where a = 4/34. 

Let us consider triangulations with two boundary loops 7 l\ and separated by a 
geodesic distance r (to be defined below). We call l\ the entrance loop and I2 the exit 
loop and it is convenient to view one of the links of l\ as being marked. For a given 
triangulation we have already defined the geodesic distance d(l, li) between a link / and 
a set of links l\ to be min//^ d(l, /'). We say that a loop 1% has a geodesic distance r to 
another loop l\ if all links / £ 1% have a geodesic distance r to l\. Note that the definition 

5 We find it convenient in this section to change from the variabie z used in the fast section to x = 1/z. 
6 The constants c t are determined by eqs. (3.29) and (3.32) of the last section. 
We use the same notation li and I2 for the set of boundary links and the number of boundary links. 
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(4.1) 
(4.2) 
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Figure 6: A typical graph contributing to G |U (/, /'; r). A blob represents a marked link. 



is not symmetric in l\ and l 2 . However, it is the natural definition for the application we 
have in mind, as will be clear from the composition law (4.8). We denote the class of such 
triangulations T(r, h,l 2 ). The geometry of a triangulation in T(r, l\, l 2 ) will be that of 
a cylinder with "height" r (see fig. 6). T(r, h,l 2 ) is related to the class of triangulations 
T(/i) with a boundary of length l\ and one marked link as shown in fig. 7. As can be 
seen from fig. 7 the elements in T(r, /1, ^) are obtained by chopping off a disk T" from the 
disk T G T(/i) such that a loop of length l 2 and geodesic distance r to /1 is created. As 
illustrated in the figure the loop l 2 might only be one of many loops in T which have a 
geodesic distance r to l\. 

We can now define a generalization of the two-point function G^r) to the two-loop 
function, where the loops are separated by a geodesic distance r: 

G,(h,l 2 ;r)= £ e"^, (4.7) 

T£T(r,l u h) 

Like the two-point function, we expect from general considerations that the two-loop 
function falls of exponentially for r — > 00 with the same exponent as the two-point function. 
In addition one has the following fundamental composition law which comes from the fact 
that a cylinder of height r\ + r 2 has a unique decomposition into cylinders of height r\ 
and r 2 , respectively (see fig. 8): 

00 

G^h,l 2 ;r 1 + r 2 ) = J2G^h,l;r 1 ) G^(l,l 2 ;r 2 ). (4.8) 
;=i 

This law expresses nothing but the summation over intermediate states since these can 
depend only on the equivalence classes of the one-dimensional (induced) metrics, and the 
equivalence classes of one-dimensional metrics are uniquely characterized by the length of 
the loop. Using (4.8) we can in principle find G |U (/i, l 2 ; r) if we know G |U (/,/ / ; 1) and it is 
seen that G |U (/, /'; 1) via (4.8) acts as a transfer matrix. It is possible to find the transfer 
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Figure 7: Illustration of how to obtain a graph in T(r, /, /') from a graph in T(l) by 
chopping off a disk with the boundary length /'. As shown there can be many loops of 
length /' and geodesic distance r from / in a specific graph belonging to T(/). 



matrix G |U (/,/ / ; 1) by purely combinatorial arguments. As for the disk amplitude and the 
multiloop correlators it is convenient to introduce the generating functional: 

oo 

G tl (x,y;r)= ]T x^y^G^h^r), G tl (x,y;r=0) = -^—. (4.9) 
h,h=i 1 Xy 

The last equation follows from the observation that 

G tl (h,l 2 ;r = 0) = 5 luh . (4.10) 
We note that G |U (/i, 1 2] r) can be reconstructed from the knowledge of G l _ t (x, y; r): 

c -<'" w - " = k ^ i ^ G '«»- " r) (4 ' n) 

where the contours C x and C y surround the origin and avoid the cuts of G l _ t (x,y;r). 
Further it should be noted that the fundamental composition law (4.8) can be rewritten 
as: 

f dz 1 

G fM (x,y;r 1 + r 2 ) = f —^G fM (x,-;r 1 )G fM (z,y;r2). (4.12) 

JC 27T1 Z Z 



This allows us to write: 



dz „ , 1 



G^x, y; r + 1) = & — - G^x, -; 1) G^z, y; r). (4.13) 

JC 27T1 Z Z 

Our main problem will be to find an explicit expression for G l _ t (x,y; 1). This problem 
was first solved in [17, 18]. 
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Figure 8: A typical graph which explains the composition law of G |U (/, /'; r). 
4.2 Solution of the combinatorial problem 

Let G^(x,y) denote the transfer matrix where the entrance loop is not marked. In the 
counting of different configurations the marking of the entrance loop will convert one 
unmarked loop into l\ marked loops, i.e. G^\x,y) is related G fl (x 1 y 1 1) by: 



The transfer matrix G^ '(x,y) is the sum of all possible graphs connecting loops separated 
by a geodesic distance r = 1. For each triangle there will be a factor g = e _A1 . In fig. 9 
we show a typical graph. One should notice that the graphs of this kind are obtained by 



(4.14) 



combining around the boundary four types of graphs, 




blobs (7) denote disk topology). In this way G^ (x,y) can be written 



as 




log(l 







(4.15) 



where the factor 1/n comes from the cyclic symmetry. The explicit expressions for //\^ 




are as follows: 




(4.16) 




(4.17) 
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Figure 9: A typical graph contributing to G^\x,y). a thin broken line represents an 
entrance loop and a thick broken line represents an exit loop. 



.A 9 Fix: q) — 1 , 

4 = gx 2 1 f , (4.18) 
Q = x 2 F(x;g), (4.19) 

since £^ contains two links from the exit loop, one from the entrance loop and one 

triangle, i.e. a factor y 2 xg and .jpj, contains one link from the exit loop, two links plus 

any disk graph with one marked vertex (where the two links join the disk amplitude) and 
one additional triangle, i.e. a factor yx 2 F(x; g)g. Similar explanations are immediately 

obtained for and . Substituting this into (4.15), we get (with g = e _A1 ) 

G { °\x, y) = - log(l - gxy 2 - gx 2 F(x;g)y - gx(F(x;g) - 1) - x 2 F(x;g)). (4.20) 

The continuum limit of the above expression is taken in the same way as when we 
discussed the continuum limit of multiloop correlators Wo(zi, . . . , z s ) in the last section, the 
variable x of the entrance loop being analogous to the variable 1/z of the loop correlators, 
as is clear from the relation between F(x; g) and Wq{z). On the other hand it follows from 
the composition law (4.12) that the continuum limit of x and y should be taken around 
values which are inverse to each other, i.e. the variable y of the exit loop is similar to z 
itself. Consequently we perform the expansion: 

x = x c (l-aX), y =—(l-aY), g = g c (l-Aa 2 ), (4.21) 

X c 

where x c = l/z c . In addition we know from the properties of the two-point function that 
the geodesic distance at the continuous level, R, should be related to the geodesic distance 
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r at the discrete level as in (2.29), i.e. introducing the lattice scaling a(pi): 

Rr^ra 2 ", (4.22) 

and that the continuum two-point function is related to the discretized two-point function 
by a multiplicative renormalization (2.32). Let us assume that we also have a multiplicative 
renormalization of the two-loop function G l _ t (x, y; r): 

G A (X,Y;R) = a s °G ll (x,y;r), (4.23) 

where Sq is a real number. As for the multi-loop correlators Wo(zi, ...,z s ) the contour inte- 
grals in the composition law (4.12) go into inverse Laplace transforms by the substitution 
(4.21) and we get: 

ioo 

/dZ 
—G A (X,-Z;R 1 )G A (Z,Y;R 2 ), (4.24) 
Ztti 

— ioo 

i.e. we find Sq = 1. Substituting (4.21) into eq. (4.20), we obtain 

G(°\x,y) = - log [const. {(X + Y)a - F A (X)a 3/2 + C(a 2 )}] , (4.25) 

where we introduced the following notation for the "scaling" part f(x;g) of the disk 
amplitude given by eq. (4.2): 

f^x) = f(x-g) = T A (Xyl 2 + 0(a 2 ) (4.26) 

T A {X) = a(X- ly^Wx + VA, (4.27) 



and where a = y2/(9\/3 — 1)- Thus, we obtain the transfer matrix at the continuous 
level, 



G A (X, Y;R = a 2 ") = aG,(x, y;r=l) = -j^y - a 1 ' 2 A 



-MX) 



X + Y 



+ 0(a). (4.28) 



Note that the leading term in (4.28) is the Laplace transformed delta function S(L — V). 
By substituting (4.28) into (4.24) for r\ = 1 and r 2 = r, we obtain 

ioo 

2 _ [ dZ I 1 u/2 d T A {X) 



G A (X,Y;R + a 2 ») = j — (^-^ - a^ 2 — -f^ + 0{a) )G A (Z, Y; R) 

— ioo 

= G A {X,Y-R)-a l l 2 -^[T A {X)G A {X,Y-R)} + 0{a). (4.29) 

We only get a non-trivial continuum limit if 2v = 1/2, i.e. the Hausdorff dimension 
d h = l/v = 4. 

If we use v = 1/4 we obtain the continuum differential equation [17]: 

— G A (X, Y;R) = -— [T A (X)G A (X, Y; R)] . (4.30) 
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Figure 10: Fundamental peeling decompositions, (a) removes a triangle while (b) and (c) 
remove a two-folded part. The solid line with the mark represents the entrance loop and 
the dashed line the exit loop after one peeling step. 




(b) 



Figure 11: Decomposition of the surface by (a) slicing and (b) peeling. 



Let us check the universality by deriving the differential equation (4.30) differently. At 
the discretized level the procedure outlined above can be viewed as a "slicing decomposi- 
tion" of the triangulations. Starting from the entrance loop we "slice" the triangulations 
in cylinders of thickness r = 1. Let us instead consider a different set of "deformations" 
which can bring us from the entrance loop to the exit loop, the so called "peeling decom- 
position" [18]. In this decomposition we start at the marked link at the entrance loop and 
remove the triangle containing this link. We might have a situation where the marked link 
is not associated with a triangle, but is folded with another link. In this case we simply 
remove both links. These three fundamental steps in the peeling decomposition are shown 
in fig. 10. Successive application of these steps is like peeling an apple (see fig. 11), in 
contrast to the situation above which was like slicing the apple and the individual steps 
in the peeling decomposition can be viewed as l/l of the slicing decomposition, where / is 
the length of the boundary. 



27 



The fundamental equation satisfied by G |U (/, /'; r) by a one-step peeling will be: 

1-2 

G^l, /'; r) fl G^(/ + 1, /'; r) + 2 £ F(Z"; </)G^(/ - /" - 2, /'; r). (4.31) 

i"=o 

If we identify the deformation (4.31) with (l//)-step of the slicing decomposition, we find 

-— G^l, /'; r) = gG^l + 1, /'; r) + 2 £ F„ (/")<?„(/ - /" - 2, /'; r) - G |U (/, /'; r). (4.32) 
' ^ r /"=o 

Taking the discrete Laplace transform of (4.32) we get: 



d_ 

dr 



G^x^-^) = x—(2x 2 f ll (x)G ll (x,y,r)j, G tl (x,y;r = 0) = ^ _ . (4.33) 



In this case it is manifest that one can obtain (4.30) after taking the continuum limit [18]. 

The peeling decomposition has a number of advantages compared to the slicing decom- 
position. One can easily use it to construct a string field theory [18] and from the point of 
view of this article it has the advantage that we have a simple differential equation even 
at the level of discretized loop length. This allows us to avoid certain ambiguities of first 
taking the discretized loop length to infinity by the substitution (4.21) and next taking it 
to zero in order to construct the two-point function. 

4.3 Solution of the differential equation 

We now discuss the solution [8] to eqs. (4.33). It is readily found since we have a first 
order partial differential equation: 



x fii{x) xy 



z 2 fn(x) 1- xy' 



(4.34) 



Here x(x,r) is the solution to the characteristic equation of the partial differential equation 
(4.33). The integral of the characteristic equation is 



-( x ' r ) dx' 



2x*f li (x') 

and this expression can be by inverted to give: 



— smh W- - 

d V 1 - c 2 x' 



x'=x(x,r) 



c 2 c 2 sinh 2 (V + sinh" 1 - 5 2 ) + 5 2 ' 

where So, S\ and S 2 are all positive and defined by 

So = 9 -^(c 2 -c 1 )(c 2 -co) = 0((A f i)i), 



Si 
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(c 2 - Cj)(c 2 - c ) 
C2(C1 - C ) 

co(c 2 - Ci) 



C2{C\ - C ) 



(4.35) 
(4.36) 

(4.37) 
(4.38) 

(4.39) 
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It is readily checked that x — > 1/ 1 ci for r — > oo and x(x, r = 0) = x. 

In principle we can calculate G^(/i,/2;r) from eqs. (4.12), (4.34) and (4.36). Here let 
us only verify that the exponential decay of G^(/i,/2;r) is independent of l\ and l^. For 
r — t- oo one gets 

G^(h, l 2 - r) = const. 8 S 1 e~ 2S ° r + C(e" 45or ). (4.40) 

where const, is a function of 0(1) which depends on Co, c\, C2,h and l^. 

We can express the two-point function G ^(r) in terms of the two-loop function and 
the one-loop function. Let us consider a marked link. For a given triangulation we can 
systematically work our way out to the links having a distance r from the marked link 
by peeling off layers of triangles having the distances 1, 2, . . . , r to the marked link. After 
r steps we have a boundary consisting of a number of disconnected boundary loops, all 
with a distance r to the marked link. One of these is the exit loop described by the two- 
loop function and we get the two-point function by closing the exit loop of length I2 by 
multiplying the two-loop function G |U (/i = l,/2;r) by the one-loop function 8 i^/2) an d 
performing the sum over 1%, i.e., as shown in fig. 12, 



G>M = E^Gi = M2; r)^^) 

Ox Jc y 2my y Oy 

1 d 



x=0 



2=0 g d 



-F, 



r 



x=0 



(4.41) 



As long as ci — c\ is small and r is larger than a few lattice spacings, we get [8]: 

cosh(^or) 
L sinh 3 ($or 



G,(r) = const.^i— - TTT^ i 1 + O(5 )) . (4.42) 



Formula (4.42) shows how to take the scaling limit: Let us return to the original formu- 
lation and write in the limit Afj, — > 0: 



cosh (Au) 4 3r\ 

G,{r) = const. (A M ) 3 / 4 ± ^-L. (4.43) 

sinlr (Afj,)if3r\ 



where const, and fi are positive constants of 0(1) {fi = y/6g c )- 
We conclude the following: 



1 



1. G>(r) falls of like e 2 {^^) 4 P r f or r _^ 00) i. e . the critical exponent v = | and the 
Hausdorff dimension d# = 4. 

2. Gfj,(r) behaves like r~ 3 for 1 <C r <C A/2~i, i.e. the scaling exponent 77 = 4. 

1 

3. The pre-exponential factor to e _2 ( A ^) 4 ^ r is (A^)^ -1 , i.e. precisely the factor needed 
in order to take a continuum limit according to the general discussion following 
eq. (2.27). 



8 To be more precise we have to multiply the two-loop function G M (7i,l2; r) by I2 since the exit loop is 
unmarked and we can glue the marked one-loop cap to it in I2 ways. 
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Figure 12: The 2-point function represented as a summation over 2-loop functions times 
1-loop functions. 



4. From Fisher's scaling relation we get 7 = z/(2 — 77) = —1/2. This well known result 
can of course also be derived directly from 

00 

X(n) = G »( r ) = const - " c 2 (A^ + • • • (4.44) 

r = l 

by use of (4.43), but it should be clear that the explicit calculation in (4.44) is 
nothing but a specific example of the general calculation used in proving Fisher's 
scaling relation. What is somewhat unusual compared to ordinary statistical systems 
is that the anomalous scaling dimension r] > 2. r] = is the ordinary free field result, 
while r] = 2 is the infinite temperature limit, and for statistical systems we expect 
i] < 2. 

5. Any two-loop function G^(/i,/2;r) has the same behavior as G^r) as long as l\^2 
stay finite as A/2 — > 0. 

It is clear that we could have taken the continuum limit almost at any point in the above 
calculations and we get 

G(R; A) = Hm^^M ~ (A) 3/4 co S h[(A)^] ^ 

sinh 3 [(A)4i?] 

where we have rescaled R with a factor (3 compared to the definition R = r a 2 " in eq. (4.22). 
The factor in front of G^r) is the usual "wave function renormalization" present in the 
path integral representation of the propagator and the unusual power M = 2A 1 / 4 of the 
"mass" is due to dfj = 4. 
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5 Euclidean quantum gravity in d > 2 



5.1 Basic questions in Euclidean quantum gravity 

The moment we address Euclidean quantum gravity in dimensions higher than two a 
frightening number of basic questions appears. Let us just list some of them: 

(1) : How do we cure the unboundedness of the Einstein-Hilbert action in d > 2? 

(2) : Does the non-renormalizability of the gravitational coupling constant not break any 

hope of making sense of the theory? 

(3) : What is the relation between Euclidean and Lorentzian signature in the absence of 

any Osterwalder-Schrader axioms to ensure that we can rotate from Euclidean space 
to Lorentzian space-time? 

(4) : What is the role of topology, keeping in mind that four-dimensional topologies cannot 

be classified? 

One could hope that a non-perturbative definition of quantum gravity could help us 
in understanding the first three points. We propose a such a definition [19, 20, 21] here. 

Not much is known about the question of summing over topologies. As long as we think 
in terms of continuum physics and write down the path integral it offers the possibility 
of summing over different manifold structures as well as integrating over inequivalent 
Riemannian structures of a given manifold: 



In two and three dimensions we do not have to worry about the meaning of Top, since there 
is equivalence between smooth manifolds and topological manifolds. In two dimensions the 
manifolds are uniquely characterized by their Euler characteristic \ an( i the summation 
over Top is simply a summation over \ or the genus g = 1 — x/2 of the surfaces. In spite of 
this simple prescription surprisingly little progress was made in defining the sum in (5.1) 
using continuum methods. The matrix model formulation of simplicial quantum gravity 
at least allowed us to address the question in a quantitative way by means of analytic 
continuations of the integral (3.8), but so far the results were ambiguous. If we move to 
three dimensions we encounter a slight classification problem in the sense that no simple 
parametrization of the various topologies exists. However, the problem gets completely 
out of control when we move to four dimensions. For four-dimensional manifolds there 
is no equivalence between smooth and topological structures. Topological manifolds exist 
which do not admit smooth structures and some topological manifolds admit infinitely 
many inequivalent smooth structures. If we insist on summing over all smooth structures 
J2top wm t" e rather unwieldy. To complicate the operational meaning of J2rop it should 
be added that four-dimensional manifolds are not algorithmic classifiable, i.e. no finite 
algorithm in the sense of Turing exists which allows us to decide if two arbitrary four 
dimensional manifolds are equivalent. On the other hand arguments (not known to us) 
might exist which dictate a restriction of the allowed class of manifolds. Since there seem 
to be fermions in the world one could argue that the manifold should be a spin manifold. If 
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one makes the additional (rather arbitrary) restriction that it should be simply connected 
such a (smooth) manifold is characterized by its Euler number and its signature where 
these in four dimensions are given by: 



X(M) 



1 



d 4 ?V^ T^abcdJLa'Vc'd' 



aba'b' ^cdc' d' 



(5.2) 



t(M) 



128tt 2 
1 



/ d 4 ^TZ abcd TZ ab c , d ,£' 
Jm 



M 



cdc'd' 



(5.3) 



967T 2 



For simply connected spin manifolds the signature r is a multiple of 16, while \ is an 
integer > 2. Thus eq. (5.3) seems a minor extension compared to two dimensions where 
the summation is over Xi but the restriction to simply connected spin manifolds is not 
natural at this stage of a quantum theory of space time. 

While these problems tend to discourage any attempt to make sense of the path integral 
of Euclidean quantum gravity, it is still our obligation to try to investigate if it is possible. 
Below we will argue that the use of dynamical triangulations, which work so well in two 
dimensions, allows us to discuss several of the above issues in some detail and might be a 
candidate for non-perturbative definition of quantum gravity even in higher dimensions. 

5.2 Definition of simplicial quantum gravity for d > 2 

Let us define simplicial quantum gravity as a generalization of the two-dimensional con- 
struction in two dimensions: In d dimensions (where d = 3 or d = 4) we construct all closed 
(abstract) simplicial manifolds from K rf-dimensional simplexes. As in two dimensions we 
imagine that the lengths of the links in the simplexes are a (which we take as 1 unless 
explicitly stated). For such a combinatorial or, equivalently, piecewise linear manifold we 
can apply Regge calculus and in this way assign a Riemannian metric to the manifold. By 
such an assignment we see that the discretization is able, for finite K , to assign a meaning 
to the sum and integral in (5.1) 



The discretization has the same virtue as in two dimensions: In principle it allows a 
unified treatment of the summation over topologies and Riemannian structures. 

Few remarks should be said about the formula (5.4). First one could try as in two 
dimensions to make the gluing automatic, and in this way arrive at generalized matrix 
models [27]. Until now these models have not been as useful as in two-dimensions, the 
reason being that although the models allow a 1/N expansion, the two-dimensional in- 
terpretation as an expansion in topology is absent. Secondly, in light of the complicated 
relation between topology and diffeomorphism for four-dimensional manifolds one could 
be worried that similar problems arise when we compare combinatorial, i.e. piecewise lin- 
ear, manifolds and smooth manifolds. However, for dimensions d < 7 we have equivalence 
between piecewise linear and smooth structures. Whatever subtleties might be involved 
in defining the sum on the lhs of eq. (5.4) it should be captured at the rhs of eq. (5.4). Of 
course eq. (5.4) itself is rather formal as it stands. The problem, which exists even in two 
dimensions, is that if no restriction is imposed on topologies the number of triangulations 




(5.4) 
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grows at least as fast as Kl, K being the number of triangles. No reasonable discretized 
action can kill this entropy factor and 

Z = J2e~ ST (5.5) 

X 

will be divergent. The naive expression (5.5) makes no sense. In two dimensions the 
way out was to fix topology. This bounds the number of triangulations exponentially and 
eq. (5.5) will be well defined. Only afterwards the topology was allowed to fluctuate 
and a summation attempted in the context of matrix integrals. In higher dimensions it 
is sensible, as a minimum, to try to define eq. (5.5) for a fixed topology. We need the 
following theorem 

Theorem: The number of combinatorially equivalent rf-dimensional manifolds is an ex- 
ponentially bounded function of the number of d dimensional simplexes. 

We call two simplicial complexes combinatorially equivalent if they have a common sub- 
division and when we talk about equivalence classes of piecewise linear manifolds we have 
in mind combinatorial equivalence. While the theorem has been known in two dimensions 
for many years and it has been rather convincingly established numerically in three di- 
mensions [22], there have recently been some controversies concerning its validity in four 
dimensions [23], but now there appeared a proof which is valid in any dimension and for 
any fixed topology [24]. 

For d > 4 manifolds exist which are not algorithmic recognizable. As a curiosity we 
can mention that for such manifolds their number is not algorithmic calculable. This does 
not mean that the number cannot be exponentially bounded, only that there is no finite 
algorithm which allows us to calculate the exact number. 

After these general remarks it is natural as a first exploratory step to fix the topology 
of our four-dimensional manifold to be the simplest possible, that of S* 4 . The Einstein- 
Hilbert action for a simplicial manifold can be calculated by Regge calculus. However, we 
do not need the full machinery in our case where the simplexes are identical and all link 
length equal. The Regge version of the Einstein action is the sum over deficit angles of 
the d — 2 dimensional sub-simplexes times their d — 2 dimensional volume. In our case 
the deficit angle associated with arf-2 dimensional sub-simplex is 2ir — c ■ 0(^-2), 
where 0(^-2) is the order of 11^-2, i-e. the number of d dimensional simplexes of which 
nd-2 is a sub-simplex, and c is a constant. It follows that 

\U d ^K\ oc ( 27r " c-o(n d _ 2 )) ■ (5.6) 

^ J > Regge nd _ 2 

If we note that the number of d — 2-dimensional sub-simplexes in a rf-dimensional simplex 
is d{d+ l)/2 we can write: 

S[A,G] = jd^^-j^K) S T [k d ,k d _ 2 ] = k d N d {T)- k d _ 2 N d . 2 (T), (5.7) 

where the piecewise linear manifold is defined by the triangulation T and N d (T) and 
Nd-2(T) denote the number of d- and d — 2-dimensional simplexes in the triangulation T. 
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We can view \jk,i- 2 as a bare gravitational coupling constant. At first sight the action 
might seem much too simple to have anything to do with gravity. Our point of view will 
be the opposite: The fact that the action is so simple reflects the beauty and simplicity of 
quantum gravity, and hopefully this simplicity will be reflected in the solution of the theory. 
Our final prescription will be: 

Z[KG] = I* Vl7*T) e ~ Si9 '' A ' G] Z l k *-*> k <A = £ e- s rt<-» k < (5.8) 

54 ^ 1 ' Tes d 

We will be interested in d = 4 and sometimes in d = 3, in which cases we get 

Z[k u k 3 ] = E e- fc3JV3(T)+fclJVl . (5.9) 
Tes 3 

Z[k 2 ,k 4 ] = e- fc4JV4(T)+fc2JV2 . (5.10) 

xes 4 

For d = 3,4 there are only two independent coupling constants as long as we only want 
an action which depends on global quantities like A^. It follows from the so called Dehn- 
Sommerville relations which express that the rf-dimensional simplicial neighborhood of 
any sub-simplex is homomorphic to a rf-simplex: 

Ni = p(-l) k+d ( k i ll)N k . (5.11) 

If we include higher derivative terms in the action we will certainly loose the simplicity of 
eqs. (5.9) and (5.10). The higher derivative terms will contain explicitly the order of the 
sub-simplexes which carry the curvature. We will not discuss the lattice implementation 
of these. 

Let us now discuss the phase diagram of the discretized theory. Assume d = 4. Since 
it is easy to prove that N 2 (T) < c ■ N^T), where c is some constant, the conjecture 
above implies that for each k 2 , a k1(k 2 ) exists such that the lhs of (5.10) for a given k 2 is 
well-defined for k^ > k1(k 2 ) and divergent for k^ < k1(k 2 ). A potential continuum limit 
should be taken as k^ — > k\ from above. The corresponding phase diagram is shown in 
fig. 13. If for a fixed k 2 approaching k\ we shall be probing the infinite volume limit of 
the discretized system. It does not imply that there necessarily will be a continuum limit. 
Rather we should view the system as a lattice system where the infinite volume limit is 
taken. For some specific values of the bare couplings critical points might exist where a 
correlation length diverges and where a continuum limit exists. Such a point is tentatively 
indicated at the figure. Approaching this point in a specific way will then define the 
renormalized cosmological constant and the renormalized gravitational constant. Since 
we are in unchartered territory one should be open-minded for other possibilities, e.g. the 
possibility that a whole range of k 2 s might correspond to topological gravity where the 
metric, and correspondingly concepts like volume and divergent distances, play no role. 

5.3 Observables 

It is possible to define the same critical exponents for higher dimensional simplicial quan- 
tum gravity as in two-dimensional quantum gravity. 
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Figure 13: A hypothetical phase diagram for four-dimensional gravity. 



First we can define an entropy or susceptibility exponent 7. In the following it will 
always be assumed that the topology is spherical, i.e. the combinatorial manifolds are 
combinatorially equivalent to the boundary of a 5-simplex. According to the theorem 
mentioned above the number N(N±) of triangulations which can be constructed from N 4 
4-simplexes is exponentially bounded. Let us now fix k 2 . As already argued there will be 
a critical point k 4 (k 2 ). This implies that 

Z(k 2 ,N 4 ) = E e k * N * ~e k ^ N *f k2 {N 4 ) (5.12) 

TeS 4 (N 4 ) 

where fk 2 (N 4 ) stands for subleading corrections. If the subleading correction is power-like 
we define 7(^2) by 

Z(k 2 , N 4 ) ~ Nl {k2) - 3 e^ N *(l + 0{l/N)). (5.13) 

However, it is not clear that we have such a behavior. It is possible to imagine an asymp- 
totic behavior like: 

Z(k 2l N 4 ) ~ e «4 C W4-c(fc 2 )JV-(fc 2 )+-( 1 + 0{l/N)) (5.14) 

where < a < 1. In this case the exponential correction given by a will always dominate 
over the power-like correction determined by j(k 2 ). There are strong indications that 
there are several regions with different asymptotic behavior, depending on the value of k 2 . 
This will be discussed below. 

Apart from the entropy- or susceptibility exponent 7 we can introduce the critical 
exponents v and r] already discussed in detail in the context of two-dimensional quantum 
gravity, where they were determined from the properties of the two-point function. Pre- 
cisely the same construction as in two dimensions can be carried out in higher dimensions 
and the definitions of Hausdorff dimension etc. is basically the same. 

5.4 Results of numerical simulations 

Presently there has not been much progress in analyzing (5.10) by analytic methods except 
for d = 2. However, the action is well suited for the use of Monte Carlo simulations and 
the results to be discussed later have been obtained by such simulations. For details about 
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Figure 14: The average radius of the universes of sizes 9000 (triangles), 16000 (squares) 
32000 (pentagons) and 64.000 (stars) versus k^. 



the most recent simulations we refer to [26]. Here let us only discuss the results. They 
can be summarized as follows [20, 21, 25, 26, 27] 

(1) : There seem to be two different regions, as a function of the bare inverse gravitational 

coupling constant k^- For small or negative values of ki the typical quantum universe 
will be very crumpled, with almost no extension and a very large, if not infinite 
Hausdorff dimension, while the universes for large values of ki will be elongated 
with a Hausdorff dimension as small as two. In fig. 14 we have shown the average 
radius for universes of size 9000, 16000, 32000 and 64000 4-simplexes as a function 
of &2- The two phases are separated by a phase transition which is of order two or 
higher. At the transition point, k^, the Hausdorff dimension might be finite (the 
precise value is not well determined, but it could be close to four). 

(2) : The same results are valid for three dimensional simplicial quantum gravity except 

that the phase transition seems to be of first order, rather than of higher order [27]. 

From fig. 14 it is seen that the change between the elongated region and the crumpled 
region becomes increasing visible as the size of the system increases. In addition the 
critical point seems to move to higher values of ki- We observe a so called pseudo critical 
point k^N^). The indication of convergence to a limiting value k\ as the volume N4 — > 00 
is shown in fig. 15 and we conclude that we have a genuine phase transition. From the 
general theory of finite size scaling the pseudo critical point k^N^) of a first order phase 
transition scales to the limiting value k\ as I/N4. Since we observe a scaling like l/y/N^ we 
conclude that the transition most likely is a higher order transition. In three dimensional 
simplicial quantum gravity it is difficult to perform this kind of measurement since we 
observe pronounce hysteresis around the transition point. This is a typical sign of a first 
order transition. 
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Figure 15: The pseudo critical point ^2(^4) 



Let us first discuss the measurements in the elongated phase, i.e. for k 2 > k\. An ideal 
quantity to measure in the computer simulations is the number of four-simplexes (n(r)) N 
within a spherical shell of thickness 1, a geodesic distance r from a marked four simplex, 
the average taken over all spherical triangulations with one marked four-simplex. It will 
depend on the coupling constant k 2 and it is related to the two-point function precisely 
as in two dimensions: 

G(r, N 4 ) ~ JV? ( * 2) - V^)"* (n(r)) N4 , (5.15) 

provided the subleading corrections are power-like. In fig. 16 we show the measured distri- 
bution ( n(r) /N4) N just above the transition to the elongated phase and a corresponding 
fit to r exp( — r 2 /N4) . We see a very good agreement and it becomes even better if we 
move further into the elongated phase. Now recall the general scaling relations derived 
in the context of two-dimensional quantum gravity, but valid in any dimension if a finite 
Hausdorff dimension exists: 

M r ))iv 4 ~ r4h ~ l for Kr<7V 4 1/<ift , (5.16) 

(n(r)) Ni ~ r a e -c(r<H/N 4 )W h -l) foj . N l/d h < r < ^ (5 _ 1?) 

we conclude that dh = 2 (and a = 1). In addition we can measure the critical exponent 7 
very conveniently by baby universe counting. Again the arguments are identical to ones 
presented in the two-dimensional case. The result is shown in fig. 17 and it is natural from 
the figure to conjecture that 7 = 1/2 as long as we are in the elongated phase. From (5.15) 
and (n(r)) N we know G(r, N4) and we can construct the two-point function Gk 4 (r; A^)) 
by Laplace transform (as in the two-dimensional case): 

00 

G k4 (r; k 2 ) = J2 N- 3 / 2 e- A ^ N * re"^^ ~ e ~' cr ^ , (5.18) 
JV 4 = 1 
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Figure 16: The measured distribution (dots, the error bars smaller than the dots), and 
the fit (the curve) using the functional form r exp( — r 2 /N^). 



where A (^4) = ^4 — ^4(^2) is assumed to be small. This function is precisely the two-point 
function of the so called branched polymers, which are known to have internal Hausdorff 
dimension dh = 2. We conclude that the numerical simulations provide convincing evi- 
dence that the elongated phase of simplicial quantum gravity corresponds to a well known 
statistical theory, the one of branched polymers. The tendency to create baby universes 
is so pronounced in this phase that the geometry degenerates to the generic lowest dimen- 
sional fractal structure possible, i.e. that of branched polymers. 

When we lower the value of ki and move below the critical point k\ the fractal structure of 
our ensemble of piecewise linear manifolds changes drastically. A glance on fig. 14 shows 
that the average radius hardly changes with the volume. This is an indication that the 
Hausdorff dimension is large or maybe infinite. If we move deep into this phase the average 
curvature is negative and in addition there are only few baby universes and they are small. 
This could lead to the idea that we entered a phase with "smooth" manifolds of negative 
curvature. For such manifolds one would expect that the volume of geodesic balls of radius 
r would grow exponentially with the radius, which is what we observe. Clearly this is a 
"fake" infinite Hausdorff dimension and indeed we should observe the dimension dh = 4 
in the sense that (n(r)) N ~ r 3 for small geodesic distances. A closer look at "typical" 
members of the computer generated manifolds indicates that they cannot be considered 
as "smooth" . Rather they have a few vertices of very high orders which connect to almost 
any other vertex in the manifold and in such a situation it is not surprising that the linear 
extension will be small. In addition we have not been able to fit to any sensible power 
dependence r dH ~ 1 for small r. A plot of log (n(r)) N shows indeed a linearly growing 
function of r up to some ro(N4) which is not much different from the average value (r) N . 
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Figure 17: The measure 7 in the elongated phase for various size lattices (N4 = 9000, 
16000, 32000 and 64000). 



A reasonable fit to {r) N is 

(r) N4 = a(k 2 ) + b(k 2 ) log N 4 . (5.19) 

This again gives some support to the idea that the Hausdorff dimension is infinite in this 
phase since it appears as a limit of (r) N ~ N l l dH for dn — > 00. Finally the extrapolation 
of (5.17) to infinite Hausdorff dimension indicates that in such a case one should observe 

(n(r)) Ni ~ Cl e~ m ^ r (5.20) 

down to quite small distances r N^ dh , — > 00, which naturally is replaced by r > 
a\ log N4 + a 2 . This is indeed what we measure. 

The observation that (n(r)) N grows exponentially from r ~ 6 out to r ~ ro and then 
falls off exponentially indicates that we deal with an infinite Hausdorff dimension at all 
distances and it is easy to get a quite good "phenomenological" fit to (n(r)) N which 
incorporates both these features by choosing e.g.: 

(n(r)) Ni ~ exp (-mi(& 2 )r - c 2 e- m ^ k2 ^ . (5.21) 
It will grow like e( C2m2_mi ) r for small distances and fall off like 

M r ))iV 4 ~ c i e~ mi ^ r - C2 e-( mi ( fc2 ) +m2 ( fc2 ) r + • • • (5.22) 

for large distances, while a N4 dependence in the coefficient c 2 would explain the observed 
N4 dependence of ro- The data and a fit of the form (5.21) are shown in fig. 18 for k 2 = 1.26 
and N4 = 64000, i.e. right below the transition to the crumpled phase, where the fit is 
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Figure 18: exponential fit (curve) of the form (5.21) to the measured n(r; N4) (dots, error 
bars less than dot-size) in the crumpled phase (N4 = 64000, ki = 1.26). 



worst. But even so close to critical point (5.21) works quite well over the whole range of 
r > 6. It should be mentioned that the coefficient in front of the second exponential in 
eq. (5.22) is negative. This implies that the term cannot be given the interpretation as 
an additional heavier mass excitation. However, just looking at the long distance tail the 
distribution (n(r)) N allows us to determine nil, and ci from (5.22). On the other 
hand we can determine c^va^ — m\ from the short distance exponential growth alone and 
find good agreement. This indicates that long and short distance behavior are intertwined 
in the case of infinite Hausdorff dimension, as they are in the case of finite Hausdorff 
dimension, where dfj appears both in the short distance and long distance expression for 



We conjecture that the internal Hausdorff dimension is infinite for ki < k\. 

In accordance which this conjecture we have in the elongated phase no "mass" term 
mik^) = (k4 — k1(k2))" which scales to zero as we approach ^4(^2)- The exponential 
coefficient m\(k2) is finite in the infinite volume limit k4 — > ^4(^2)- However, it is most 
interesting that m\(k2) scales to zero as ki — > &2j the phase transition point between the 
crumpled and the elongated phase. This gives a strong indication that the system at the 
transition might have a finite Hausdorff dimension, which could very well be larger that 
the generic dh = 2 found in the elongated phase. In fig. 19 we have shown the scaling of 
the mass m\(k2) as a function of ki- 

The above mentioned numerical "experiments" 9 suggest the following scenario: The typ- 
ical quantum universe, determined without any Einstein action (i.e. ki = 0) has (almost) 
no extension. Its Hausdorff dimension might be infinite and internal distances between 

9 See [26] for a detailed discussion of the most recent results. 



(n(r)) Ni (see e.g. (5.16)-(5.17)). 
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Figure 19: The behavior of the masses mi^) and mi^) + ^2(^2) from (5.21) in the 
crumpled phase as a function of k\ — k<i- 



points are always "at the Planck scale" . By that we simply mean that no consistent scaling 
can be found which will be compatible with finite continuum volume and finite Hausdorff 
dimension. For a finite value of the bare gravitational coupling constant there is a phase 
transition (second or higher order) to a phase with a completely different geometry with 
pronounced fractal structures. It is tempting to view the transition between the two kind 
of geometries as a transition where excitations related to the conformal mode are liberated, 
since large ki is a region which formally corresponds to small values of the gravitational 
coupling constant. Right at the transition it seems as if we have the chance to encounter 
genuine extended structures with a finite Hausdorff dimension. Maybe the fact that the 
transition between the two types of geometry is of second (or higher) order can be used 
as the starting point for a non-perturbative definition of quantum gravity. 

Of course it is crucial to be able to perform high statistics simulations at the critical 
point in order to investigate this possibility in greater detail. 

6 Discussion 

We have shown how it is possible to discretize reparametrization invariant theories and 
apply with success the methods known from theory of critical phenomena. In two- 
dimensional gravity it was possible to solve the theory and calculate the Hartle-Hawking 
wave functional and reparametrization invariant observables like the volume-volume cor- 
relator as a function of the geodesic distance. These observables are not easily calculable 
using a continuum framework like Liouville theory, not to mention canonical quantization. 

In three- and four-dimensional gravity we also obtained a well defined theory with 
specific scaling relations and it was possible to determine these reliably by numerical 
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methods in part of the parameter space, the parameters being the bare gravitational 
and the bare cosmological coupling constants. An interesting phase transition point was 
found in four dimensions and it is a potential candidate for the point in coupling constant 
space where a non-perturbative definition of quantum gravity might lead to a physically 
acceptable theory. 

It is interesting to compare the results obtained with theoretical suggestions. The 
simplest comparison is made with so called "conformal gravity". Presently there are two 
versions which are quite different: The first one is based on adding higher derivative terms 
to the Lagrangian and taking the coupling constant in front of the square of the Weil 
tensor to infinity [28], while the other approach appeals to the scale anomaly as the factor 
which dictates the important physics [29]. Both theories have a non-trivial fixed point. In 
the first case it is an ultraviolet fixed point, in the second case an infrared fixed point, so 
physics is quite different for the two theories. They both give definite scaling predictions 
for the partition function and it seems as if there is surprisingly good agreement between 
the predictions coming from the theory with an infrared fixed point and simplicial gravity 
at k\. However, the comparison is not unambiguous and it is too early to make definite 
statements. This work is in progress. 

At the moment we let loose topology we also loose control of the theory. In two 
dimensions the use of dynamical triangulations led to a number of very interesting attempts 
to perform the summation over topologies by means of analytic continuation of the matrix 
integrals which can be used to represent the partition function. It is a most important 
task to understand if it is also possible in higher dimensions to address in a quantitative 
way the summation over topologies. Dynamical triangulations will hopefully help us to 
address the question and force us to ask the right questions, since the formalism is so 
simple, that it is impossible to hide in mathematical abstractions. 
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